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Abstract 

We point out that the inflaton inevitably couples to all non-conformally coupled 
matters gravitationally through an oscillation in the Hubble parameter or the cosmic 
scale factor. It leads to particle production during the inflaton oscillation regime, 
which is most efficient just after inflation. Moreover, the analysis is extended to 
the model with non-minimal inflaton couplings to gravity, in which the Hubble 
parameter oscillates more violently. We apply our results to the graviton production 
by the inflaton: gravitons are also produced just after inflation, but the non-minimal 
coupling does not induce inflaton decay into the graviton pair. 


1 Introduction 


Inflation explains both the homogeneous and isotropic universe and the small density 
fluctuation in a beautiful way. In a successful inflation model, the energy density of 
the inflaton must be converted into high-temperature radiation. This process is called 
reheating. 

In a standard picture, the reheating is caused by the perturbative decay of inflaton. If 
we introduce a coupling between inflaton and light matter, which may or may not be the 
standard model particles, the inflaton has a hnite lifetime and it decays at some epoch 
in the early universe. Since we do not know detailed properties of the inflaton, we can 
arbitrarily choose the inflaton-matter couplings as long as they do not disturb the inflaton 
dynamics during inflation. 

In this paper we show that inflaton necessarily couples to all scalar helds gravitationally 
through the cosmic expansion. We explicitly show that, during the inflaton oscillation 
regime, the Hubble parameter and the cosmic scale factor oscillate with (twice) a frequency 
of inflaton. That is, the scale factor depends on 0^, which leads to an effective coupling 
like C ~ 0^((9x)^/M|, where 0 is the inflaton, y is any scalar held and Mp is the Planck 
scale. This leads to “decay” or “annihilation” of the oscillating inhaton into y particles. A 
similar phenomenon, called gravitational particle production, was investigated in Ref. [T] . 
While Ref. [T] studied particle production in the transition from the de Sitter space to 
the radiation or matter dominated universe, we point out that particle production always 
occurs in the inhaton oscillation regime. Although this gravitational interaction is not 
strong enough to complete the reheating, still it may be a subject of phenomenological 
interests. Interestingly, gravitons are also necessarily produced by this process, although 
its abundance will be below the observable level of future gravitational wave detectors. 

We further discuss the case of non-minimal inhaton coupling models: f{4>)R gravity. 
In such models, the Hubble parameter and the scale factor oscillate more violently with 
time during the inhaton oscillation regime. We will explicitly show that the scale factor 
linearly depends on 0 and hence there arises a gravitationally induced coupling like C ~ 
(j){dxy/Mp. This interaction itself can complete the reheating. The same result was 
obtained in Ref. [2] in a diherent approach. Our method makes it clear that inhaton 
oscillation explicitly couples to the matter sector inevitably. 

As a by-product, we will show that the inhaton 0 does not always decay into y 
particles even if it has a coupling like 0(0y)^/M in the Lagrangian for f{4>)R models. 
Depending on the choice of M, the 0 decay into y particles may be forbidden due to 
some non-trivial cancellation. Actually, this cancellation necessarily happens in the case 
of inhaton-graviton coupling. 

Before going into the main text, we comment on the terminology of “gravitational 
particle production” often used in diherent contexts and their applicability. 

• Originally Ref. [1] called the particle production due to the non-adiabatic change of 
the background at the end of inhation as the “gravitational particle production”. 
This is close to what we are going to study, although our main focus is the inhaton- 
oscillation dominated era after inhation. We will show that this may be interpreted 
as the “decay” or “annihilation” of the inhaton through the gravitational interaction. 
In our viewpoint, it is clear that particles lighter than the inhaton are necessarily 
produced. 
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• Quantum fluctuations generated during inflation for nearly massless scalar fields 
will eventually become a coherent oscillation mode, which is also called as “gravita¬ 
tional particle production” [1]. This production mechanism is applicable for scalar 
helds whose mass lighter than the Hubble parameter during inflation and we do not 
consider this production mechanism in this paper. 

This paper is organized as follows. In Sec. [2] we study the system of Einstein gravity 
and minimally-coupled matter and show the gravitational coupling between inflaton and 
ordinary matter helds. In Sec. [3] we extend our analysis to the case of non-minimally 
coupled inhaton. Sec. 0] is devoted to summary and discussion. 


2 Einstein Gravity 


In this section we discuss particle production after inhation in the standard Einstein 
gravity. Let us consider the following action 






--Mr - 


+ ^i 


M, 


( 2 . 1 ) 


with the Robertson-Walker metric, 


= —dt^ -I- a^{t)dx^, 


( 2 . 2 ) 


where g = det (g^j,^), a is the scale factor, R is the Ricci scalar, cj) denotes the inhaton, Sm 
collectively represents any other held and we assume that S'm does not explicitly contain 
the 0 held. We will see how the inhaton 0 couples to helds in S'm in the inhaton-oscillation 
dominated regime in this minimal setup. 


2.1 Background dynamics 

Let us assume that the inhaton dominates the universe. We consider the homogeneous 
background dynamics of the inhaton. The Friedmann equation is 

3MM = (2.3) 

where H = d/a is the Hubble parameter and = 0^/2 -|- V. The equation of motion of 
0 is 

0 + 3iL0 -t-100 = 0, (2-4) 

where = dV/dcf). From the equation of motion, we obtain 

00 + 3iL02 _ 0. (2.5) 


Here, we denote 0 as a deviation from the potential minimum. 

In the deep oscillation regime m 0 S> id, where m 0 = is the inhaton mass 


scalej^ P 0 is a good conserved quantity in an inhaton oscillation time scale. Therefore, 


from the Virial theorem, we obtain 


(02)=n(R), 

Note that this condition corresponds to 0/Mp <C 1. 


( 2 . 6 ) 
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where we have assumed power-law potential V (x. (jf' and (• ■ ■) means time average in a 
period longer than oscillation period but much shorter than the Hubble time scale. Then 
we obtain 




(2.7) 


This shows that (p<^) oc and {H) = {n + 2)/(3nt). 

This oscillation-averaged picture is often useful when discussing cosmological evolution 
of the oscillating inflaton held. However, it does not contain small oscillations in and 
H, which come from the effects neglected in the Virial theorem. Although the oscillation 
amplitude is small, it induces particle production and may have observational impacts in 
some cases® 

To analyze the oscillatory behavior, quantities are divided into oscillation-averaged 
parts and oscillation parts: P(j, = {p^) + dp^j, and H = {H) -|- 5H. The oscillation part 
satishes 


5H 


3n 

n + 2 


{Hf 


2Mj,' 


( 2 . 8 ) 


By noting 0^/2 + V= (p^) -|- 6p^ ~ 3Mp (77)^ (1 -|- 26H/(H)), this is rewritten as 


• 677 

6H+ - {H)6H ~ 

n + 2^ ' 


1 1 
n + 2M], 



(2.9) 


The hrst term on the l.h.s. is 6H ~ m^6H 3> {H)SH, and hence it dominates the second 
term. For the same reason, on the r.h.s., the hrst term gives the dominant contribution. 
Then, by integrating this equation, we obtain 


6H ~ 


n + 2 Ml' 


( 2 . 10 ) 


The scale factor a{t) is obtained by integrating a/a = H = (H) + 6H: 


a{t) ~ (a(t)) 1 - 


1 02 _ 


2(n + 2) 


Ml 


{a{t)) -ai[j 


n + 2 
3n 


( 2 , 11 ) 


where quantities with subscript i are evaluated at an arbitrary initial time in the inhaton- 
oscillation dominated era® This explicitly shows that the Hubble parameter H as well as 
the scale factor a{t) depend on 0^ and oscillate with time. Note that the relative oscillation 
amplitude, 5H/H cx 0/Mp, becomes smaller and smaller as the universe expands since 
the inhaton oscillation amplitude decreases with time. Thus the oscillation in the scale 
factor is most efficient just after inhation. 

We have performed numerical calculation to check the above formula. Fig. [T] shows 
time evolution of H (left) and 6H (right) for n = 2. We have compared numerical results 
and the approximate analytic formula fl2.10p for 6H. We have taken 0j = 0.1 and = 1 
in Planck unit. 


*2 Such kind of particle production is also briefly discussed in ref. [3]. 

Here we neglect terms suppressed by ^ {(fP)/Ml in the oscillation-averaged scale factor (a(t)). 
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Figure 1: Time evolution of H (left) and 6H (right). We have compared numerical 
results and the approximate analytic formula fl2.injl for 6H. We have taken = 1 in 
Planck unit. 


2.2 Gravitational annihilation of inflaton 

2.2.1 Scalar 


Now let us discuss the gravitational inflaton decay. To be concrete, we consider the matter 
Lagrangian for a real scalar y: 


Sm= 


-^g^'^df.xduX - ^rnlx^ 


( 2 . 12 ) 


where we assume m^. As shown above, the scale factor and hence contains 

dependence. Therefore, neglecting terms including m^, the action can be expanded as 

Sm = J drd^x {a{t)f (^1 - \ “ i^ixf] , (2.13) 


where we have used the conformal time dr = dt/a{t) and the prime denotes derivative with 
respect to r. This explicitly shows that the inflaton 0 couples to {dxY and 0 (partially) 
“decays” or “annihilates” into y particles. According to the analysis of particle production 
under the oscillating background 0 M, it might be interpreted as the annihilation of the 
inflaton into y particles. Thus we call this “gravitational annihilation” for convenience in 
the following. 

To estimate the production rate, we write down the equation of motion for Xk = CLXk 
where Xk denotes the Fourier mode of y with comoving wavenumber fc. 



Xk — 0 


Xk + 


- 2 {nf + 



Xk = 0 , 


where V, = a' / a - 
creation with k ~ 


= aH. Therefore, the effective mass of y oscillates rapidly@ 
rritj) occurs and the creation rate was studied e.g. in Ref. [5]: 


r(# -t yy) ~ 


C 171^ 
3277 Mj, Mj, ’ 


(2.14) 

Particle 


(2.15) 


If y couples to the Ricci scalar as £ = ^f{x)R, the equation of motion becomes y" — V^y + 
{a"/a){x + = 0- Thus if /(y) = —y^/6, y does not couple to the scale factor or the inflaton in the 

limit rriy^ = 0. This corresponds to the conformal coupling. 
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for n = 2 where C is a numerical constant of 0{1), with $ being the inflaton oscilla¬ 
tion amplitude. Since <h is decreasing, it never exceeds the Hubble parameter and the 


reheating is not completed by this gravitational annihilation J*d S till, however, this small 
annihilation process may yield signihcant amount of particles 
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Here we comment on the possibility of parametric resonance. The width of the hrst 

The Hubble 


♦7 


resonance band for at /c ~ rritj, is Ak ~ qrricj, with q = 
expansion removes the produced x particles from the instability band in the phase space 
within the time interval ~ g/if. Since the time scale for the exponential growth of the 
X particle is ~ we need 3> if in order for the parametric resonance to 

occur EE]. This condition is rewritten as ($/Mp)^ 3> 1, which cannot be satished as 
long as <h -C Mp. Thus we neglect the effect of parametric resonance and simply use the 
production rate fl2.15p . 

Now let us estimate the abundance of x particles produced by this gravitational an¬ 
nihilation. Because of the <h dependence, the abundance is dominated by those created 
just after inflation: 


n^(iii„f) ~ 


P<p r(00 ^ xx) 




H, 


inf 


9C 

1677 


H: 


inf 5 


(2.16) 


where iiinf is the Hubble parameter just after inflation. Note that this is of the same order 
as the number density found in Ref. [T] where the particle creation at the transition from 
the de Sitter universe to the matter or radiation dominated universe was studied. Our 
approach shows that particle creation always occurs in the inflaton-oscillation dominated 
era and makes it clear and explicit that the gravitational particle production is induced 
by the inflaton-scalar coupling in the Lagrangian. 

If X is massive (but much lighter than the inflaton) and stable, and assuming that its 
interaction is so weak that they cannot be thermalized, its energy density per entropy 
density is given by 


^ ^ 9C m^TuHinf 
s 6477 Mp 


~ 8 X 10 -®GeVC 


( '^x \ 

U 06 GeV/ 


T, 


R 


IQIO GeV 


/ f^inf A 
VlQi^GeVy ’ 


(2.17) 


where Tr denotes the reheating temperature, which is completed by some other inflaton 
interaction terms. If x is absolutely stable, this can be one of the production mecha¬ 
nisms of present dark matter. On the other hand, if x is a moduli which decays into 
standard model particles through Planck-suppressed interactions, this is another source 
of the cosmological moduli problem 

For nearly massless x particles such as axion-like particles, the gravitationally pro¬ 
duced X particles may contribute to dark radiation. The energy density of such massless 

It is pointed out that, if reheating ends only through the gravitational particle production in the 
context of ref. [I] , the graviton overproduction can be problematic [7] . 

Note that $ in Eq. (12.151) should be regarded as amplitude measured from the potential minimum 
of the inflaton. 

The width of the instability band may be interpreted as follows. If one regards the particle production 
process as an annihilation (pcf) —>■ xx> the dispersion relation of x is given by = rn^{l —(jp/2Mp). Since 
b is oscillating, the momentum of x has a finite width of Ak = m0d)^/2M|> = qm^. 

In the context of gravitational particle creation of Ref. [T], see also Refs. 00 ■ Ref. [5] studied the 
moduli abundance in the form of coherent oscillation induced by the non-adiabatic change of the Hubble 
parameter. 
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X particles, normalized by the radiation energy density, is estimated as 


3C o-iy r V''7 

Vprad/T=r 167 r Mj, VlQi^GeV/ VlOi°GeVy VlO^^GeV 

(2.18) 

where is the Hubble parameter at the completion of the reheating. Thus its contri¬ 
bution to the present dark radiation energy density is negligibly small. 

One comment is in order. A scalar particle often obtains a Hubble-induced mass 
through the coupling with Ricci scalar R or the Planck-suppressed interactions with in- 
flaton held. If such a mass term is absent, i.e. if y is nearly massless during inhation, 
it necessarily develops quantum huctuations of overhorizon modes which in turn will be¬ 
come a coherent oscillation after inhation, and this may be a dominant contribution to 
the X abundance [1]. On the other hand, if x obtains a Hubble mass during inhation, its 
quantum huctuation is suppressed. However, even in such a case, our production mecha¬ 
nism always works as long as S> iPinf, which in fact many inhation models including 
new inhation and hybrid inhation satisfy, and this may be a dominant contribution to 
the X abundance. In this sense, the estimate fl2.17p gives a robust lower bound on the x 
abundance unless x is conformally coupled to the Ricci scalar. 



2.2.2 Vector boson and fermion 

Next let us consider gravitational particle production of massless vector bosons and 
fermions. In order to see this, we hrst consider the following action for a massless vector 
boson: 


Sm = -\J (2.19) 

where F^^ is the held strength of the vector boson. If we use the conformal time dr = 
dt/a(t), the action can be rewritten as 

Sm = -\J dTd?xg^^g'^^F^,F^0, (2.20) 

where is the metric in the Minkowski space-time and the summations are taken over 
{r, Xi, X 2 , 0 : 3 }. From this, it is obvious that the massless vector boson does not couple 
to the scale factor. Note that this discussion holds also for non-abelian gauge helds. 
The same thing happens for the case of a massless fermion. To see this we consider the 
following action: 


Sm = - j d'^xe'ipe'^'j’^Vf.'ijj, (2.21) 

where is the spinor held, is the vierbein and e = det(e“) = y/-^- Here we use 
the Greek letters p, z/,... as the Einstein indices and the Latin letters a,b,... as the local 
Lorentz indices. Moreover, is the covariant derivative which is dehned as 

v, = d, + ( 2 . 22 ) 
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where the spin connection is given as 

^ (2.23) 

By using the conformal time and rescaling the spinor held as '0 = the action can 

be rewritten as 

Sm = - j (2.24) 


where the summation for fi is taken over {r,xi, 2 : 2 ,xs}. Therefore, the scale factor does 
not couple to the massless fermion. Note that these facts are closely related to the Weyl 
invariance of massless vector bosons and fermions. In other words, this annihilation 
process is applied to any helds lighter than the inhaton if they are not Weyl invariant 
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For example, if vector bosons and fermions are massive, they are coupled to the scale 
factor and hence produced by the inhaton oscillation. 


2.2.3 Graviton 


Finally, let us see the inhaton-graviton coupling. The graviton is dehned as the transverse- 
traceless part of the metric perturbation: 

ds^ = —dt^ -I- a{tY{5ij -I- hij)dx'^dx\ (2.25) 


where hij satishes dihij = ha = 0. Expressing the two helicity modes by (A = -1-, x), it 
is known that the equation of motion of the graviton is the same as that of the minimally- 
coupled massless scalar held. Thus hx = ahx satishes the same equation as 02.141) . It 
means that the graviton production rate is given by 


T{(j)(() —)■ hh) ~ 


C 

IQti Ml Ml' 


(2.26) 


The abundance of graviton, or the gravitational wave, is given by Eq. 02.181) . The present 
frequency of the gravitational wave is 


/gw — 


m0a(Ffinf) 

27r 


~ 2 X 10® Hz 




1013 GeV/ VlOioGeV 


T, 


R 


1/3 


H, 


inf 


1014 GeV 


-2/3 


(2.27) 


where we have taken the present scale factor to be equal to one. The present gravitational 
wave spectrum in terms of flew has a peak at this frequency and scales as for higher 
frequencies. Around this frequency range, the abundance of gravitational wave is typically 
below the observable level of future space laser interferometers [II]- 


3 Extended Gravity 

In the previous section we have shown that, in the Einstein gravity, the inhaton necessarily 
couples to scalar helds and graviton through the 0 dependence of the scale factor, which 
is proportional to In this section, we discuss particle production in the extended 

See Refs. [TU] for the gravitino production. 
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gravity, since it sometimes happens that the Hubble parameter as well as the scale factor 
oscillates more violently [T2l[T3]. This means that the scale factor linearly depends on 0. 
Therefore, we will expand H and a(t) at the linear order in 0 and neglect 0^ dependence 
in this section. 

Let us consider the following action: 




J 


- 7:9^'"dt,(t)d^(t) - H(0) 




M, 


(3.1) 


where /(0) is an arbitrary function of 0, which can be expanded as 


/(0) — Mp + ci-^ H-(3.2) 

Here we take /(O) = Mp without loss of generality by redehning 0 as a deviation from 
the potential minimum. 


3.1 Background dynamics 

Let us hrst consider the homogeneous mode of 0 which undergoes a coherent oscillation 
around the potential minimum: 0 = 0(f). Let us compute 0-dependence of the scale 
factor a{t). From the Einstein equation, we obtain 

3fH^ = p^- 3Hf, (3.3) 

where = 0^/2 -|- H and 

{3H^ + 2H)f + f + 2Hf 

The equation of motion of 0 reads 

^ + 3H<P + V^- 31^{2H^ + H) = 0. (3.5) 


1 V 


+ H. 


(3.4) 


One of these three equations is actually redundant. 

It should be noticed that the Hubble parameter H is an oscillating function with 
time, as opposed to the case of the Einstein gravity. To see this, we expand H and 0 
as H = Hq + Hi and 0 = 0o + 0i, where the subscript 0 represents the solution in the 
Einstein gravity limi t f( 0) = Mp, and we regard the effect of non-minimal coupling ci as 
a small perturbation l^^u Thus we have Hq = (n + 2)/{3nt) for V oc (jH. We also consider 
the deep oscillation regime: i.e., H <C with being the mass scale of 0. 

From Eq. fl3.3p . we obtain 


Hl = 

Hi = 


P4>o 


3M|,’ 

P4>i 

QMlHo 


Cl 

2Mp 




(3.6) 

(3.7) 


where p^i = 0o0i + In order to estimate the magnitude of p^i, we rewrite fl3.5p as 


p(j>i + 6iLo0o0i + 3iLi0Q — 3ciMp^o[2Hq H) — 0. 


Note that this expansion is different from what we have done in the previous section. 


(3.8) 











This is further rewritten as 


P(j)i + 6ifo0o0i + 




2ci 


P <1)0^0 ■ 


( 3 . 9 ) 


The second and third terms on the l.h.s. are of the order of ~ and are subdominant 

compared to the hrst term. Thus integrating this equation, we obtain c^MpHl4>o, 

since 0 is an oscillating function with time scale ~ 1/m. Therefore, on the r.h.s. of 
Eq. fl3.7p . the second term dominates over the other terms. As a result, we obtain 

H, ^ ( 3 . 10 ) 


Integrating a/a = H, we hnd the scale factor 


a{t) ~ ai 




ao(t) 1 


Cl 0 A 

2 Adp J 


( 3 . 11 ) 


These results show that the Hubble parameter H as well as the scale factor a explic¬ 
itly depend on 0(f) and hence rapidly oscillating functions with time. The combination 
a?{t)f{(/)) becomes 

a\t)f{</)c^Mla,{tf. (3.12) 

Thus the 0 dependence cancels out in the combination a‘^{t)f{(j)). 

We have performed numerical calculation to confirm the above considerations. Fig. [2] 
shows time evolution of 0(f) (top left), H (top right), (bottom left) and a?/{(j)) (bottom 
right) for n = 2. We have taken 0j = 0.1, ci = 0.3 and rntp = 1 in Planck unit. We have 
compared numerical results and approximate analytic formula fl3.10p for H and fl3.1ip for 
a^. It is clearly seen that H as well as the scale factor a oscillates with time, but the 
combination a^/(0) does not. 


3.2 Gravitational decay of inflaton 

3.2.1 Scalar 


First let us consider the inflaton decay into a real scalar y. We introduce the following 
term in the action 



-^h{<t>)9‘“'dicXdyX 


( 3 . 13 ) 


where h(0) is an arbitrary function of 0, which is expanded as h{(j)) = 1 -|- di(j)/Mp 
The reason for the introduction of this non-minimal coupling will become clear soon. 
Naively, this term leads to the decay rate of the inflaton as r(0 —)■ xx) = (1287rMp). 

But this naive expectation is not correct as we will see below. 

Clearly h(0) is an oscillating function and it might seem that it induces particle cre¬ 
ation of y. But we must be also careful on the prefactor x/—g, which could cancel the 
oscillation of h(0). More precisely, by using the conformal time and metric matter action 
becomes 

Sm = J dTd^xa‘^{t)h{(j))^ [y'^ - {dixf] . (3.14) 
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Figure 2: Time evolution of 0(t) (top left), H (top right), (bottom left) and a^f{(j)) 
(bottom right). We have compared numerical results and approximate analytic formula 
fId.lOp for H. We have taken ci = 0.3 and = 1 in Planck unit. 
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Thus it is the combination a‘^{t)h{(j)) that determines the particle creation rate of x- By 
using the result obtained in the previous section, the combination a^(t)h(0) becomes 


a^{t)h{(j)) ~ ao{ty 


1 + (di — Cl) 


Mr 


(3.15) 


Thus the 0 dependence vanishes for di = ci, i.e., /(0) = Mph(0). In this case, 0 does 
not couple to x, hence no x particle production is induced by the 0 oscillation except for 
the effect studied in the previous section. In the present toy model, there is no reason to 
expect /(0) = Mph(0) and hence the decay occurs unless there is a tuning, but later we 
will see that this tuning automatically happens if we regard x as the graviton. 

Having revealed the 0-x coupling, we can calculate the decay rate of 0 as 


r(0 ^ xx) = 


ml 


1287rM2 


(di - ci)^. 


(3.16) 


It coincides with the result in Ref. [2] for di = 0. It is noticed that even for di = 0, i.e., 
no explicit coupling between 0 and x, the rapid oscillation of the scale factor leads to the 
efficient decay of the inflaton. 


3.2.2 Vector boson and fermion 


Next let us consider vector bosons and fermions: 


^M = 


d‘^Xy/^ 




(3.17) 


where hc{(t)) and hp{(j)) are arbitrary functions of 0, which can be expanded as hc{4>) = 
1 + dc^/Mp + • • • and hp{(j)) = 1 + dpcfy/Mp + • • •. As shown in Sec. 12.2.21 by using the 
conformal time and appropriate rescaling of the helds, they can be rewritten as the form 
that does not include the scale factor: 


^M = 


drd^x 




(3.18) 


Thus the oscillation in the scale factor does not affect the decay rate of 0 into massless 
vector bosons and fermions. The decay rate into vector bosons is given by 

r(0 ^ gg) = -^^dlml, (3.19) 

where N denotes the number of vector bosons. The decay rate into massless fermions 
vanishes. 


3.2.3 Graviton 

Here we apply our results to the inflaton-graviton coupling. Assuming that the back¬ 
ground is homoge neo us and expanding only the tensor part as Eq. fl2.25p . we obtain the 
graviton action 

S = J dTd^xa{tff{(j))^ [hfj - {dihijf] , (3.20) 

Terms proportional to (hij)^ is cancelled out by using the background equation (13.41) . 
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where hij satisfies hu = dihij = 0. As shown in Sec. m the combination a{tyf{(f)) does 
not depend on <p. The 0-dependence of scale factor a(t) through the Friedmann equation 
cancels the 0-dependence of /(0). Hence there does not appear an infiaton-graviton 
coupling like 0(0h)^ in /(0)i? models. 

3.3 Coupling of the fluctuation 

In Sec. 13.2.11 we have seen that the homogeneous coherently oscillating scalar couples to 
scalar field x unless /(0) = Mph(0). In this section we consider the same matter action 
as that in Sec. 13.2.11 Let us assume that 0 sits at the potential minimum 0 and consider 
the fluctuation 6(j)(t,x) around it. Hence the background part 0 satisfies the following 
equations: 


3LfV(0) = ^(0), (3.21) 

VS) - = 0. (3.22) 

If V (0) = 0, we are dealing with Minkowski space, while H(0) > 0 corresponds to de Sitter 
space. We want to know the coefficient of coupling like ~ 50(0x)^ in the Lagrangian. 
This is a non-trivial task because of the mixing between 50 and the fluctuation of scalar 
part of the metric due to the non-minimal coupling. 

We make use of the ADM formali sm M to deal with the metric fluctuation and 
specify the physical degrees of freedom l^^ 

= —N‘^dt^ + '^ij idx'' + I3^dt) {dx^ -f- jd^dt) , (3.23) 

where N is the lapse function, 0* is the shift vector and is the metric of space. The 
inverse metric is given as 


-1 0 * 

0* AT^yi - 0*0t 

In this parametrization, the Ricci scalar can be decomposed as 

R = fl® + - E^) + 2V^ f - 0AP>]V. (3.25) 

where the superscript (3) denotes the covariant quantity with respect to the spacial metric 
7 jj. In addition, Eij is related to the extrinsic curvature as 


(3.24) 



Eij = NKij = 



(3.26) 


and E = Y^Eij. Moreover, is the unit normal vector of the time-like hyper-surface: 

n" = i (1, -/3‘) . (3.27) 


The analysis of fluctuations during inflationary regime using the ADM formalism is done e.g. by 
Ref. [H]. 
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By using Eq. fl3.25jl . the Lagrangian can be rewritten as follows: 


^ 


f(4.) i (£«By - B") 


E 


-/'(0)^ (0 - 0*0,0] - iV7*^Vf 0,/(0) 

^ ^0 - 0 * 0 , 0 ) ' - y 7*''0*00,0 - NVicfy) 


2N 




2N 


(3.28) 


The equations of motion for the lapse function and shift vector give the constraints among 
the helds. The constraint equations derived from the Lagrangian fl3.28p are 


lapse : 0 =^m («<=» - - E^)\ + - f)'di4>) 




2N^ 


shift : 0 =V 


(3) 


N 


fmE]-6)E) 


0 - 0 * 0 * 0 j --7*^ 0 . 00 ,- 0 -^( 0 ), 

+ -^djcf) 


N 


7 ( 3 ) 


N 


(0 - 0 * 0 . 0 ) 


N 


(0 - 0*0i0)0,-0, 


(3.29) 


(3.30) 


where we omit ^-dependent terms by assuming that y is a small quantity. 


413 


We now expand the scalar held and metric around the background part as follows 


414 


N — It o, 

0* = 0T, 0.0X = 0, 

7ij = e^‘^+^^(e'*)y, dihij = ha = 0, 

0 = 0 T 50, (3.31) 

and treat a, 0^, Ci 50 as perturbations. The background scale factor is given by 

a{t) = e'*. It is enough to solve the constraint equations fl3.29p and fl3.30p to hrst order 
in the perturbations in order to obtain cubic terms in the perturbations. This is because 
contributions from higher order terms vanish due to the background equations of motion. 
Then, for our purpose, it is enough to expand /(0) to hrst order in the perturbations: 

/(0) - /(0) + = Ml(^ + , (3.32) 


where we set /(0) = Mp and /,</,(0) = c^Mp in the second line. Only terms up to this 
order is relevant for the constraint equations. It is straightforward to solve Eqs. 03.291) 
and 03.301) by substituting the perturbative expansion 03.311) . The solutions are 


a = C 


Cl 50 



(3.33) 


413 7 pjg treatment is justified since we solve the constraint equations only to first order in perturbations. 
It completely fixes the gauge condition for the perturbations. 
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where we used the background equations of motion 03.221) to obtain this result. This 
relates the metric fluctuations a and ( to 6(j). 

Then turning back to the Lagrangian fl3.28p . the last term is expanded up to hrst 
order in the perturbations as 


S= d‘^xa{t)^- 


<fi-^ + 3C 


Mf 


d^xa{tf^{di- Cl) ^ 


X 


a ) X - 

(da? 




(3.34) 


This agrees with the (p-x coupling obtained in the previous subsection in the case of 
homogenous oscillating 0 held. In particular, the coupling vanishes for di = ci. Moreover, 
we can see that the coupling of the scalar huctuation with gravitons also vanishes. In 
the ADM parametrization, the graviton action in the form of scalar-graviton-graviton is 
written as 


S = f d'^x 


Cl 


6(p 


+ 3C 


a 


(hij) 


Cl 


6(p 


+ C + Ct 


{dih 


i]) 


(3.35) 


By substituting the solution of the constraint equation fl3.33p . we hnd that the scalar- 
graviton-graviton interaction vanishes. Therefore, the non-minimal coupling like f{4>)R 
does not induce the inhaton decay into gravitons. This result is consistent with that of 
the previous section. However, we note that there still remains a coupling 02 (ah)2 as 
studied in Sec. [21 


3.4 Conformal transformation 


The calculations so far indicate that the 0 decay into x P^h is forbidden if /(0) = 
Mph{(j)). These results may be understood as a conformal transformation. We have 
analyzed the toy model in the Jordan frame, but we can move to the Einstein frame in 
which the gravitational action takes the Einstein-Hilbert form by the following conformal 
transformation: 


Qeflu = (3.36) 

Here the subscript J and E denote the quantities in the Jordan and Einstein frame, 
respectively. In this subsection, we keep the subscript J to make the discussion clear. By 
the conformal transformation, the Ricci scalar transforms as 

Rj = D2 + QOp In D - In In D). (3.37) 

Then the action fl3.1l) combined with fl3.13p is rewritten as 


S= d'^xyj-g-p, 


-MpRp — -gp — Ve(0e) — x 


1 Mph(0j) 


2 /(0j) 




where 


Re(0e) — 


Vjjh) 


d0E = Mf 


'2/+ 3/2^ 


2/2 


d0j. 


(3.38) 

(3.39) 

(3.40) 


Thus it is evident that the 0-x interaction vanishes for /(0) = Mp/i(0). 
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4 Summary and Discussion 

We have shown that the inflaton necessarily conples to any scalar held as long as it is not 
conformally conpled to the Ricci scalar throngh the cosmic expansion during the inhaton- 
oscillation regime. The point is that the Hubble parameter and the cosmic scale factor 
have small oscillatory features with (twice) the frequency of the inhaton oscillation. It 
induces an effective coupling like C ~ 0^((9x)^/Mp that leads to the “annihilation” of 
the inhaton coherent oscillation into x particles. Although the rate of particle creation is 
rapidly decreasing with time compared with the Hubble parameter and hence the reheat¬ 
ing is not completed by this gravitational annihilation, still it can have phenomenological 
consequences such as the heavy dark matter production or the cosmological moduli prob¬ 
lem. Here we stress that any scalar particle lighter than the inhaton is produced by this 
process, even if it has a Hubble-induced mass or a mass larger than the Hubble scale dur¬ 
ing inhation. The same particle production also applies to the graviton; high frequency 
gravitational waves are necessarily produced after inhation. 

Moreover, we have investigated the case of inhaton non-minimal coupling: f{4>)R 
model. In this model, the Hubble parameter and the cosmic scale factor oscillate more 
violently with the frequency of inhaton oscillation. That induces an ehective coupling like 
C ~ (j){dxy /Mp that leads to the decay of inhaton into y particles. The decay rate can 
be large enough to complete the reheating. However, we also found that such a linear 
coupling to the graviton is automatically cancelled. 

Our hnding may be understood as particle production in a time-dependent background 
in a broad sense. In the context of preheating EE] it is well known that the explicit 
coupling like C ~ or leads to efficient particle production. We emphasize that 
even if there are no explicit couplings among light particles and the inhaton, they actually 
couple to each other gravitationally. More formally, the factor y/—g = a{t)^ in the action 
is directly related to the inhaton dynamics through the Friedmann equation, which may 
also be regarded as the constraint equation that connects the metric to the inhaton. In 
the Einstein gravity, the oscillatory behavior of the scale factor is rather mild, but in some 
extended gravity models, the oscillation might be more drastic so that the gravity-induced 
couplings are strong enough to complete the reheating. 

Although in this paper we have considered the inhaton, similar discussions are applied 
to any coherently-oscillating scalar held, because the scale factor oscillates even if the 
energy density of the scalar held is subdominant in the universe. 
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A Generalization of the Inflaton Sector 


In Sec. 13.31 we consider the couplings between the fluctuation 6(j){t,x) and x for the 
minimal kinetic term of 0. The same analysis can be done for a more general form of 0 
sector. In this section, we consider the following Lagrangian: 







h(0) 

2 


g^'^df.xduX , 


(A.l) 


where X = Note that this Lagrangian includes a variety of extended 

gravity models such as f{R) models. 

As in Sec. 13.31 let us assume that 0 sits at the potential minimum and consider the 
fluctuation 6(j){t,x) around it. The background part 0 satisfies the following equations: 

0 = K^{^,0) + 6H^lS)^ (A.2) 

0 = K{^,0) + 3H^f{^). (A.3) 


By using the ADM parametrization, the Lagrangian can be rewritten as 


S = / 


f(4,) i - £=) 


{i> - P'di<l) - iV7«vf>8j/(,« + NK(,I,,X) 


2N 2 

The constraint equations derived from the Lagrangian flA.4l) are 

lapse : 0 A® - 

K, 


- + K{4,, X)-‘-^[4,- P'dt 


shift : 0 =V 


(3) 




+ 


UE 


3 


^(0 - 0 * 0 ^ 0 ) 


N 
Kx 
N 


a0 


(0 - 0*0i0)0j0, 


(A.4) 


(A.5) 


(A.6) 


where again we neglect y-dependent terms here. We now expand the scalar held and 
metric around the background part as follows: 


X — 1 T cn, 

0* = 0^, W = 0, 

Xij = dihij = ha = 0, 

0 = 0 + 50, (A.7) 

and treat a, 0^, h, 50 as perturbations. The background scale factor is given by 
a{t) = e^. Here we solve the constraint equations flA.bli and flA.bp to first order in the 
perturbations. In order to do this, notice that the iL^-dependent terms can be neglected 
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since they are second order in the perturbations. Moreover, we can take K{(I),X) = 
K{(j),0) = V{(1)) since X is a second order quantity. We again expand /(0) to hrst order 
in the perturbations: 


/(0) - /(0) + /,<a( 0)^0 = M|, 1 + Cl 


50 

Wp 


(A.8) 


where we set /(0) = Mp and /,</,(0) = ciMp in the second line. By noting these facts, it 
is straightforward to solve Eqs. flA.Sp and flA.bp . The solutions are 


a = C = 


Cl 50 


2 Mp 


/^T ~ 0) 


(A.9) 


where we used the background equations of motion (IA.2p and flA.Sp . Note that these 
solutions are the same as those in Eq. fl3.33p . Then turning back to the Lagrangian flA.4p . 
the last term is expanded up to hrst order in the perturbations as 


S= d‘^xa{ty 


410+ 3C 


Mp 


= / d xa{t)'^-{di — ci) 


5(t) 


a I X - 

(day 


4^ + C + a 


- 


(A.IO) 


This coupling is totally the same as that of the minimal case. In particular, the coupling 
vanishes for dp = cp. 


B Calculation in Another Gauge Condition 

For the check of our calculation in Sec. 13.31 in this section we take a different gauge 
condition from Eq. (13.3ip . We expand the helds as 


Af — 1 T cr, 

0* = + 0T) ^i/^T = 0, 

Xij ® 0 ) 

0 = 0 + 50, (B.ll) 


and treat l3p,hij,6(j) as perturbations. In the same way as before, the constraint 
equations fl3.29p and (13.301) can be solved to obtain 


a = 


0 = 


Cl 


d 


2MpH dt 




Cl 


,-2cr 


2MpH 


50, 



(B.12) 
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We consider (j){dh)‘^ couplings in this gauge. The relevant part can be written as 

^ ciM|. 6(j) / 

8 Mp y 

Ml ■ \ 

^ ^ ^ijhij T e dihijdihijj 

- -^di'ipdi'yij'^ij. (B.13) 

By substituting the solutions flB.12p and doing some integration by parts, we dually obtain 
the following result: 


ij hij 


e'^dihiMik 


Hj 




where 5Chh/ 
as 




C\ Scj) ’ SUfifi 

I i"i ■i 


5<j)hh 


2HMp 5h,i 




(B.14) 


1st 


is the variation of the kinetic term of the graviton Chh-, which is given 


Ml■ • 

Rhh ^ ^ hijhij e dihj^jdik 




(B,16) 


The term remaining in Eq. flB.ldp can be absorbed by the following held redehnition: 

Cl 6(j) 


kij hij 


2H 


(B.16) 


After such a held redehnition, there is no (j){dhY. Therefore, a scalar held cannot decay 
into gravitons in this gauge. This is consistent with our previous result. In fact, the 
solutions in the present gauge flB.12D and flB.1611 are related to those in the previous 


gauge 03.331) by a gauge transformation 
calculation. 


4kl5 


This completes the consistency check of our 
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